NASA TECHNICAL MEMORANDUM 


NASA TM- 77720 


STABILITY OF FLUID FLOW IN A CYLINDRICAL ANNULUS 


Hubert Ludwieg 


liA.'A-Tfl-7? )) STABILITY OJF 1L0ID JPLCi IB B85-10300 

A CI.INDBICA- ABBQLOS (National Aeronautics 
and S;ace Adninistratxon) 17 p 

HC L U/&F A01 CSC! 20D Dnclas 

63/34 23303 


Translation of "Stabilitat der Stromung in einem 
zy lindrischen RingrauirT , Zeitschrift fur Flugwissenschaften, 
Vol. 8, No. 5, May I960, pp. 135-140. 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON, D.C. 20546 AUGUST 1984 



ST AMO A VD TlTCf rtoi 



X fUclplM*# €•«#(#• 

i TtiltiMbMrit 

STABILITY OF FLUID FLOW IN A CYLINDRICAL 
ANNULUS 

9. (Mt 

Auciust 19 84 

9. fifltmlai Or|«filNrt«i C «*• 

t. AyAHl) 

Hubert Ludwieg 

9. Of|«ai|«tiM Nt, 

1#. VwfcUMtN*. 

•. 

II. Ci.i»m «t 0»**» Wa 
NASw ^562 

9. PhI»*hi| Nan* m>4 AMMM 

SC IT RAN 

Box 5456 1 

111 

Tranalatloa 

12. AMNft Nom 

1 tional Aa.oaautica «o4 Spec* Adalaio tratloa 
Vaanlngtoa, L'.C. 4056# 


M. fc—t—twf Af— CM# 

19* hffl—mtpy N#m 

Translation of "Stabilitat der Stromung in einem zylindri- 
schen Ringraum", Zeitschrift fiir Flugwissenschaften, Vol. 8, 
No. 5, May 1360, pp. 135-140. 


| It. MlMtl' 

The present paper deals with the investigation of the helical flow in an 
annulus between two coaxial cylinders a regards its stability against the 
formation of helical vortices of the type known as Taylor's annular 
vortices. Assuming the annulus to be small and the velocities to vary 
linearly with radius, it is shoi.-n that the probler" can be reduced to the 
classical case of flow between two rotating cylinders. .An appropriate 
stability criterion for helical flows is derived from Rayleigh's stability 
criterion applicable to such flows. 
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STABILITY OF FLUID FLOW IN A CYLINDRICAL ANNULUS 


Hubert Ludwieg 
1. INTRODUCTION 


Consider a fully developed rotationally symmetric f lev.’ in a /135 
cylindrical annulus. In such a flow, there is no radial component 
of the velocity, and the longitudinal and tangential velocity 
components do not vary longitudinally or tangentially but depend 
only on the radius. Thus, only two types of instabilities are 
important . 

The first type of instability arises due to separation of the 
flow on the inner boundary of the cylindrical annulus and gives 
rise to the so-called dead water. This instability has been 
discussed by several authors [1,2] and will not be considered 
here . 

The second type of instability is fairly well known for the 
special case of a fully developed rotationally symmetric flow 
in an annulus m which the longitudinal velocity component 
vanishes. This type of flow occurs in the annulus between 
rotating coaxial cylinders. The case of inviscid flow has 
already been discussed by Lord Rayleigh [3] . It was shown that 
the circulation is unstable when the tangential velocity decreases 
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faster than 1/r. The instability of similar flow when 
friction is considered has been treated by G. I. Taylor [4] where 
it has been shown that friction stabilizes the flow, i.e., the 
stability limit is extended into the previously unstable fluid 
flow region. For large Reynolds numbers (where the effect of 
friction goes to zero) , Rayleigh stability limit is reached again. 

In this paper, the stability of a fully developed rotationally 
symmetric annular fluid flow is studied for the general case of 
non-zero longitudinal velocity component. Since we shall apply it 
to a fluid flow with large Reynolds numbers, we shall consider only 
the s.mple case of inviscid flow. 

2. STABILITY CRITERION 

Let us consider an inviscid fluid flow in a cylindrical 
annulus between two coaxial cylinders which is fully developed and 
rotatio.ally symmetric. The circulation is assumed to be result- 
ing from the fact that a fluid whose axial and tangential com- 
ponents depend only on r flows from a fixture V into the 
annulus R (Figure 1) . A radial equilibrium is considered to 
exist in this flow at 


v * 




Figure 1. Cylindrical Annulus 


egress (centrifugal force = force due to radial pressure gradient) . 
The fluid medium thus flows into the annulus without a change in 
velocity. We shall investigate this flow for its stability. 

(If the flow at the time of discharge from V is not in radial 
equilibrium, it is reached in a short distance from V) . 
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We shall use cylindrical coordinates r,<t>,z and denote 
the corresponding velocity components by V r , V^, Thus, 

'' 0. V, \\ (r), \\ rr- V'. (r) . (1) 

Because the flow is considered to be frictionless (inviscid) , 

V and V need not be zero on the walTs r=r. and r = r , 
q z i a 

For brevity, we set 


r i + r a = 


and denote the velocities and V 2 at the point r - r^ as 

V - and V n . For simplicity, we further assume that 

$0 zO 

r*-Ti = .ir<r t and that V, and V are linear in r. We can, then, 

$ z 

write that 
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| v, =V t , + c T (r-rJ. (2) 

• v, = V* + c, (r - rj . 


In spite of these constraints, the fluid flow described by 
Equation (2) is quite general. The flow is not irrotational and 
is not governed by constant pressure. 

Let us next consider the simple case of V z = 0. In this 
case, the flow is only in the tangential direction and so the 
Rayleigh stability criterion is valid. This criterion states 

that the fluid flow is stable when V, decreases slower than 1/r 

<t> 

or perhaps increases, for increasing radius. 

Thus, in our notation, the flow is stable for 


For, 
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perturbations known as Taylor ring vortices appear and lead 
to a quick change in velocity profile. 

The occurrence of the Rayleigh stability limit pro- 
portional to 1/r can be easily understood when a ring shaped 
fluid element coaxial with the annulus is considered. Due to 
rotational symmetry, parts of this fluid element are dispersed 
radially, so that a ring shaped element having a larger radius 
is formed. From the law of conservation of angular momentum, 
this displacement results in a change in tangential velocity 

which is proportional to 1/r. If the velocity of the outer 
fluid in the neighboring region decreases faster than 1/r, the 
centrifugal force occurring in the ring dominates the restoring 
force due to pressure gradients of the surrounding fluid flow 
and thus causes instability. If in the outward direction 
decreases slower than 1/r or even increases, the force resulting 
from pressure gradients exceeds the centrifugal force, and the 
flow is stable. 

The case ^ z =v z q » i-e., c z = 0, does not lead to anything 
new compared to the results of V z =0. We need only to cons’ ier 
the flow phenomenon in a reference frame moving with a velocity 
V z g. It is clear immediately, that we have the same flow as in 
the case of V z ^=0 , and therefore, have the same stability criterion. 

We now consider the case c 4 = 0 . We could use, here, the 
usual technique to derive a stability criterion, i.e., apply a 
suitable perturbation to the flow and with the help of the 
theory of small oscillations study the conditions under which the 
perturbation increases with time. Relative to Rayleigh case, a 
complication arises that, in addition to the generalization of 
the flow, it. is no longer justified to apply purely ring vortices 
as perturbation (Taylor vortices) . We need to also consider 
perturbation due to helical vortices. These calculations can. 
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however, be performed economically if we could, by making a good 
approximation, possibly reduce the problem to the Rayleigh case. 


We consider the flow in a reference frame moving with a 
constant velocity U along the z-axis. In this reference system, 
we have once again a fully developed rotationally symmetric flow 
in a ring annulus, having velocities and velocity-gradients 

~ ~ U, c„ c, 

respectively. The velocity U is chosen such that the slope of 
all helically shaped streamlines of this flow is independent of 
the radius r. The planes containing the streamlines which are 
normal to the cylindrical surface are then purely helical. The 
condition for U is thus given by 


V,* ~U _ V m - U + c, dr 

_ * Zji r* — — • — — — — — • 


By approximating the solution to be linear, we obtain for U/V T , 
and v^V,. 


U l^n 

1 ~c,rjv„- 

V»_ c,rJV„ 
V f . 1 - c T rJV„ 


(4a) 

(4b) 


where <5 is the angular slope of the streamlines relative to 
the radius r^ in the moving frame. 

Consider now a tube of streamlines which is formed from two 
neighboring helical planes with slope <5 and the walls of the cylinder. 
It can be seen that by changing the reference system, an almost 
planar flow is obtained. This tube of streamlines is different 
from those from the Rayleigh case only by the fact this tube is 
also curved in a direction perpendicular to the streamlines. And, 
as we advance along the direction of flow, this tube of constant 
cross-section slowly rotates around its axis. It should be realized 
that the stability or instability of the Rayleigh flow arises due 
to the interplay between centrifugal force and the pressure gradient. 

0Ri3<r:.-.- ' ■ 
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It follows immediately that, under the assumption of narrow annulus 
(Ar/ro << l) , only the curvature along the streamlines is important for 
stability and that both deviations from the Rayleigh case discussed 
above have no influence on the stability. We can, therefore, apply 
the Rayleigh stability criterion, and introduce the following 
effective quantities 


v «Wv„* + vv = c ^. 


hi = c, cos d + c, sin d . 


By setting these quantities in Equation (3), the following condi- 
tion is obtained 

c, rt r n . « _ r* (c, + c t tg A) _ ^ 

V..« ~ V'.„ 


Using the value of tangent 6 from Equation (4b) , we obtain as 
stability criterion of our flow*: 

C T r * _ ( c » V^»l)* ^ . 

1 ~c t rjV„ ( 5) 


3. DISCUSSION (OF THE RESULT) 


For reduction of the variables, we substitute for V^ 0 /r 0 in 
our stability calculation, the dimensionless velocity gradients 

r = Ct r * r _ c i r t 

’ V,.' 1 V„ 

The stability condition, then, becomes 


c - c * ^>-1 

ORIGINAL P-\r T r ’ 1 -c, 

0F POOR QUALITY 


Subsequently, it can also be verified numerically from the 
differential equations for a friction free flow that helical 
vortices which correspond to Taylor ring vortices, having velocity 
given by Equation (4a) and slope given by Equation (4b) can 
occur in the fluid flow when condition (5) is not satisfied. This 
will be shown in the appendix. 
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The resulting stable and unstable regions are drawn in a o- 
plane in Figure 2. We see immediately that for c : ~o , 
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Figure 2. Stable and unstable regions in 
? . ; plane. 

7 * 


the Rayleigh condition ?., > i is obtained again. For C 

z 

values different from zero, the condition c r > - 1 (weaker 
than potential vortex) is insufficient for stability. Thus, an 
additional z-component of the velocity gradient C functions 
always as a destabilizer, irrespective of its sign. The stability 
limit for C z values for which c, >-i, is given by the straight 
line c* = l , is shown in Figure 2. In this case, the neutral 
equilibrium is reached by a distribution of the tangential 
velocity component V^, which corresponds to a rigid body 
rotation. (V^ is proportional to r.) 

Let us consider the form of the unstable perturbation. It is 
clear immediately by differentiation that it follows the 
direction of velocity flow in the reference system moving with 
velocity U. We thus have a stationary helically shaped vortex 
whose slope is given by Equation (4b) . In a c r ,c»- diagram (Figure 
3) , the direction 6f a helical vortex corresponding to the point P 
having value c^c,- and consistent with Equation (4b) is given by 
a straight line through the points P and R (coordinates 1,0), if the 
<J> direction is considered as c<j> -axis and the z-axis as • c r . In 
a stationary reference system, these helical vortices move with 
a velocity U along the z-axis. 
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Figure 3. Change in c t . values for an unstable flow 
shown at point P. 

Let us now consider an unstable flow which flows from a 
fixture V into an apparatus, as shown in Figure 1. Then, the 
perturbations will grow very rapidly and produce a momentum 
exchange which will give rise to an effective velocity gradient 
. The components of the gradient vector (c ? ,r,) along the 
direction of the vortex decrease whereas those perpendicular to 
it remain unchanged. The latter do not correspond to shear 
but only to a rotation of the fluid. In c,, c, -diagram (Figure 3), 
the motion is therefore along the straight line through points 
P and R. 

At the stability circle (point Q) , the momentum exchange 
becomes stationary, and the velocity profile is now stable. 
However, by some combination of normal turbulence, particularly 
for large Reynold numbers, the velocity profile may possibly 
change. In addition, because of the frictional effect of the 
walls which cannot be switched off, the velocity profile will 
change slowly due to presence of the walls, even after reaching 
a stable point. 

4 . FLOW STABILITY IN WIDE ANNULUS 


The calculations presented here are valid under the assumption 
that the annulus is small (Arlr t < l) , and that the velocity com- 
ponents can be written as linear functions of r. When these 
conditions do not hold, our stability criterion is no longer 
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directly applicable. In this case, we can, however, divide the 
cylindrical annulus into a number of small ones of differential 
width dr and apply our stability criterion to each of them. 

The question, then, arises if it is possible to make a statement 
concerning the total flow. It can be safely stated that our 
stability condition is definitely a necessary condition for each 
radius of the annulus. If this is not the case for an annulus of 
given radius, then the flow in that region, as observed from a 
suitable reference system moving along the axial direction, is 
approximately planar, and unstable as shown for the Rayleigh case. 
It may, also, be assumed that this condition is sufficient, since 
for all perturbations involving a vortex of small radial extension, 
the flow is stable when it is satisfied. 


For perturbation due to vortices of large radial width, there 
is a displacement in each of the small cylinders. This, in turn, 
releases restoring forces as seen for the Rayleigh case, so that it 
may be assumed that the complete flow is stable with respect to 
these vortices. Formally stated, it implies that the flow is 
stable, if and only if, the condition 


dV v r 
dr V„ 


is valid for every radius r. 


dm;/ 


i- 


dV v r 
dr V T 




1 . 


5. POSSIBLE EXPERIMENTAL VERIFICATION 


(7) 
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In order to verify the stability criterion experimentally, 
it can be so arranged that stable and unstable flows are 
generated in an apparatus of Figure 1 and then the development of 
flow profile along z-direction is measured. For a stable flow, 
the velocity profile must remain stationary until disturbed by a 
growth of the boundary layer on the walls. On the other hand, 
a quick immediate change in velocity profile occurs over the 
entire annulus for an unstable flow. Eventually, a stable final 
position, as shown in Figure 3, is also determined from the fact 
that after reaching this distribution, there is no change or only 
a slow change in flow profile. 
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We could now think of directly observing the appearance of 
helical vortices, and thereby, verify the stability criterion. This 
would, however, fail for the apparatus shown in Figure 1, because 
after the desired flow is generated through the fixture V, a number 
of unstable disturbances occur immediately. These perturbations 
grow concurrently and, therefore, do not lead to the development of 
a distinct helical vortex. 

The stability criterion can yet be tested in a completely 
different way. We can use two sufficiently long coaxial cylinders 
which rotate independently of each other, and can, in addition, 
move relative to each other along the axial direction. In the 
middle part of the aperture, a fluid of suitable viscosity is intro- 
duced such that no turbulence occurs. The desired flow pattern 
having several values of «v?.- is obtained by a suitable choice of 
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Figure 4. Proposed apparatus for generating helically 
shaped ring vortices 

tangential and displacement velocities. The Reynolds number is 
chosen higt enough that it has very little influence on the flow 
properties. The remaining effect of the viscosity on the flow 
profile between rotating cylinders can be easily estimated by 
Taylor computation [4] because these results can be carried over 
in our case in the same manner as for the Rayleigh case. For the 
flow shown in Figure 4, the vortices should be observable just 
like the Taylor vortices between two rotating cylinders. By 
calculating the slope of the helical vortices and their displace- 
ment velocity U, we can confirm the validity of our considerations. 
Furthermore, the appearance of unstable fluid flow should also 
be recognizable by a quick change in velocity profile. 
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APPENDIX 
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The Euler differential equations for the flow in cylindrical 
coordinates are given by 


dw, 


dw. 

, dw r 


dm. 


dt 

+ w, 

dr 

+ w <p r\ ~ 

r dtp 

4^ 

r 

dz 

ss — ' 

dw. 

dw m 

dw? 

W r W<p t 

dw v 


dt 

+U? r 

dr 

±W<p - 4- 
rdtp 

*+u>, 

r 

dz 


dw j 


dw t 

dw z 

dw z 

1 


dt 

4- Wr 

9r 

+ W<p A 4 
r dtp 

■ W 'dz = 

0 

dz ' 


l d P ' 

0 dr ' 

1 dp' 

Q rr> < 


dw, w r dw v dw t 

dr r r dip dz 


Superposition of the basic flow given by Equation (1) with a 
perturbation having velocity components^, yields 

w r = v,(r,q>,z, t), 

= V». (r) + v T (r,<p,z,t), 

= V x (r) + V,{r.<f,z,t). 
p' =» P (r) + p(r,<p,z, t). 


By inserting these components in the above differential equations 
and taking into account that the basic fluid flow satisfies these 
differential equations, and ignoring all quadratic terms in 


Vrr v,,v s , 


we obtaiu 


dt 


+ V v 


d v, 
r d<p 


dVr 

r V * +V 'di 


1 dp 
g dr ' 



dV 9 

. . dt?™ 

V. 


dr 

t>r + V* *4- 
r cty 

*V,rV t 




_ . du x 

dt 

+ dr 



dv f 

, v r ( 

dvy dv t 


dr 

4- - -f 
r 

r dtp + dz 

0. 


dvy _ l dp 

dz o rdrp 

1 dp 
£ dz ' 


( 8 ) 


We now make a trial solution for the pertubation velocities, 
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V'-V'Me*"'-*-*". 

v t =v t (r)e ,, “ z ”- SH , 
P ~ p (r) e' ,a ‘ 


( 9 ) 


Here a and y are real whereas 3 can be complex. Thus, we apply 
a vortex of the type given by Taylor vortices as perturbation; 
however, it moves with a constant velocity and is helically shaped. 
Substituting Equation (9) in Equation (8) yields, 


( - () i t- V,, ‘ i + V t tx ij v, - 


2 V, . 1 dp 

7. v- — — 


ff dr ' 

I - fi i + V 9 r i + V, a i j v v + 

idV, V 9 \. 1 y.. 

+ |~ - + _ , Pt 
\ fir r f or 

( - fii +■ V T ^ i + V z <t i j v t +■ 

dV t . 1 

+ dr *'=- q *'P> 

**?' + *'+ y ii T +*ib t =0. 

dr r r v 


( 10 ) 


Let us first consider the Rayleigh case, i.e., a fluid flow without 
the 2-component and a stationary ring-shaped vortex as perturba- 
tion. We can then set 

dv 

V, =0, 1 = 0, y - 0 , .‘Hf'h =*£•» =0 


Further we can, under the assumption that the ring aperture is 
small (.lr/r,<*i), ignore o r /r relative to other terms in the continuity 
equation and approximate 


VV- VW r 


'o* 


dV v 

dr 


c< 


dVt -c 
*' dr Cz 
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Thus, we obtain 


. . 2 V .o* _ __ 1 d P 

•»' «*■ 


fii^T + 



1 . * 

h* •* 

— dtp , 

dp. 

act u j = 0 * 

dr ^ 



®wa« wt p. .... 

Of pQrto ^ t§ 

^ Quality 


By eliminating other quantities, the following differential 
equation is obtained for v r 


&b r 

dr* 



2 V, n 

fir r » 



+ 



(ID 


This differential equation yields the well-known sinusoidal 
solutions which vanish at the inner and outer cylinder, only 
when the second term in the angular bracket is negative. It 
follows, then, that the stability condition is given by 


which is identical with the Rayleigh stability criterion. 

We now consider the general case with non-zero V . By 
using 

V f = v ,t, +C <p (r ~r„), 

Vz = v. n + c t (r — r„) 


the expression -/*» + V t (,v/r)» + v,ni appearing in (10) can be rewritten 
as 

-fit + V„ y r i + V'<zi = (-#» + * +v m *)i+ 

+ ? +* c *)( r ~ r o)» +fa- 

N r O r 0 r 0 ' 

Setting 


fa + ^gx) ^ + Vm * ==* 0 7 
r 0 


(12) 
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corresponds to defining the displacement velocity fa/a and the 
slope of the helical vortex and leads to the same results as in 
Equations (4a) and (4b). Thus, we obtain 

- $i+V' r r i + V t *i- fa. (1 


By eliminating from (10) the quantities and again, omitting 

small higher order terms, and using (12) and (13) , we obtain the 
differential equation for v r 


<Pv, .|f' cA 1 , 

dr 3 * I. (f r -v; w v a )*J 

+ 2 v ,v r » ! f .*x c *_| v **)lj^® r = o. 

A*;*- • \r n n\dr 


(14) 


The stability criterion follows from (14) in a manner analogous 
to the Rayleigh case 


r„c„ - V ia/ r o ‘- Ct * ^ r7) ] “ 

= - 1 ,(Cx*+ C r * -1) >0, 


(15) 


which is identical to the stability condition (6) . 
7. SUMMARY 


Stability of an inviscid, helically shaped fluid flow in a 
cylindrical annulus is studied with respect to the formation of 
helically shaped vortices using the well-known method of Taylor 
ring vortices. Under the assumption of a small annular space, and 
tangential and axial components V and V z varying linearly with 
r, stability condition for the hetically shaped flow is derived 
from the corresponding Rayleigh case for the flow along the 
tangential direction only. Denoting the average velocity along 
the tangential direction by , the average radius of curvature 

of the annular space by r^, and the dimensionless velocity gradients 
by c, = dV, dr and o = dV,,dr with c* and ?. , the stability condi- 
tion may be expressed as 
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Hie stable and unstable regions given by this relation are 
drawn in a 2f plane in Figure 2. The Rayleigh stability 

<P Z 

limit holds for f . = o and corresponds to a tangential velocity 

component proportional to 1/r (potential vortex) . A finite 
value for c, always degrades the stability. For f t >i • the 
stability limit is given by ~* = i which corresponds to the tan- 
gential velocity component being proportional to r (rigid 
rotation) . 

The stability criterion is also given for a wide annulus. 

The validity of this criterion is, however, only made plausible, 
and has not been strictly derived. 
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